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The reaction 2 H(e, e'p)n has been studied with full kinematic coverage for photon virtuality 
1.75 < Q 2 < 5.5 GeV 2 . Comparisons of experimental data with theory indicate that for very low 
values of neutron recoil momentum (p„ < 100 MeV/c) the neutron is primarily a spectator and 
the reaction can be described by the plane-wave impulse approximation. For 100 < p n < 750 
MeV/c proton- neutron rescattering dominates the cross section, while A production followed by the 
NA — ► NN transition is the primary contribution at higher momenta. 

PACS numbers: 25.10.+S, 25.30.Fj 



For high virtuality of the exchanged photon, the 
2 H(e, e'p)n reaction is one of the simplest and best ways 
to investigate high-momentum components of the deu- 
terium wave function, possible modifications to the inter- 
nal structure of bound nucleons, and the nature of short- 
range nucleon correlations. To date this reaction was 
studied only for low Q 2 (< 1 GeV 2 ) at Saclay, NIKHEF, 
Mainz and Bates. A survey prior to 1990 can be found 
in Ref. [l[ . In general, the interpretation of these results 
suffered from large corrections due to final-state interac- 
tions (FSIs), meson exchange currents (MECs) and the 
intermediate A contribution. 

The Continuous Electron Beam Accelerator Facility 
(CEBAF) at Jefferson Laboratory (JLab) has opened a 
new frontier in the study of 2 H(e, e'p)n and other (e, e'p) 
reactions for Q 2 up to 6 GeV 2 . The first study of the 
exclusive 2 H(e,e'p)n reaction at JLab has been carried 
out in Hall A [2j. The cross section was measured as 
a function of recoil momentum, p n , up to 550 MeV/c 
in perpendicular kinematics with Q 2 = 0.68 GeV 2 . At 
low (< 300 MeV/c) recoil momentum, these data can be 
described to within l-2a by the Plane Wave Impulse Ap- 
proximation (PWIA) , while at higher momenta FSIs and 
the A contribution must be included. Two new exper- 
iments have been carried out at higher Q 2 : the first in 
Hall A for Q 2 < 3.5 GeV 2 and the second in Hall B Q 



for 1.75 < Q 2 < 5.5 GeV 2 which is reported in this letter. 
We have done a comprehensive study of the 2 H(e, e'p)n 
exclusive reaction with full kinematic coverage, which al- 
lows us to identify the dominant mechanisms. 

The experiment has been performed using the CE- 
BAF Large- Acceptance Spectrometer (CLAS) which 
consists of six sectors, each functioning as an indepen- 
dent magnetic spectrometer. Each sector is instrumented 
with multi-wire drift chambers, time-of-flight scintillator 
counters covering polar angles 8° < < 143°, gas- filled 
threshold Cherenkov counters (CCs) and lead-scintillator 
sandwich-type electromagnetic calorimeters (ECs) cover- 
ing 8° < 6 < 45°. The CLAS was triggered on scattered 
electrons identified by a coincidence between EC and CC 
signals in a given sector. 

A 5.761 GeV electron beam impinged on a target cell of 
liquid deuterium about 5 cm long and 0.7 cm in diameter, 
positioned on the beam axis close to the center of CLAS. 
The target entrance and exit windows were 15 /im Al 
foils. A 4 cm vertex cut for the scattered electron selected 
events from the central part of the target and eliminated 
events from the windows. The CLAS vertex resolution [B[ 
of a = 2 mm allowed us to estimate a background from 
the windows of < 0.5% @. 

Electrons and protons from the reaction 2 H(e, e'p)n 
were selected in fiducial regions of CLAS, where the parti- 
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cle detection efficiency is high and nearly constant. Both 
the CC and EC were used to distinguish electrons from 
pions for momenta < 2.8 GeV/c, whereas only the EC 
was used for momenta > 2.8 GeV/c where the CC be- 
came sensitive to pions. Ref. |6[ reports that n~ con- 
tamination is < 2% depending on Q 2 . The data were 
corrected for this effect. The protons were identified us- 
ing tracking and time of flight [5J] . 

The electron detection efficiency depends on the drift- 
chamber inefficiency (2.5%) and the 7r _ rejection cuts 
in the EC (2.5%) and the CC (10%), on average. The 
proton detection efficiency depends on the ir + rejection 
cut (2.5%) and the inefficiency of the drift chambers plus 
the time-of-flight scintillators (10%) [6|. 

The exclusive 2 H(e, e'p)n events were extracted from 
the data by requiring the missing mass to be that of the 
undetected recoil neutron. We measured the differential 
2 H(e,e'p)n cross section as a function of Q 2 , p n and 9 n 
(the neutron polar angle with respect to the momentum 
transfer direction), integrated it over cf> n (the azimuthal 
angle of the recoil neutron), and corrected it for accep- 
tance and radiative effects. The acceptance corrections 
were calculated using a Monte Carlo technique for all 
Q 2 ,p n and 9 n bins, and were applied event by event to 
every bin. The radiative correction factors were calcu- 
lated using the method described in Ref. Q. 

The measured cross sections (points) are shown versus 
p n in Figs. [TJ [2 and versus 9 n in Figs. [3l [4] for Q 2 = 
2, 3, 4 and 5 GeV 2 . Statistical errors only are shown. 
Systematic uncertainties due to the pion contamination, 
electron and proton detection efficiency, beam intensity 
measurements, and target density are less than 1%. More 
important are the uncertainties from the effective tar- 
get length (3.5%), acceptance corrections (5.5% point to 
point), background subtractions from the missing mass 
distributions (2-3% average; 5.5% point to point), and 
radiative corrections (4%). The total experimental sys- 
tematic uncertainty is 10% Q. 

We have investigated the same reaction theoretically 
using the most recent predictions of Ref. [8( , which have 
been programmed into a Monte-Carlo code that gener- 
ates events in the fiducial acceptance of CLAS. We sam- 
pled p n , On, 4> n , 4>e (the azimuthal angle of the scattered 
electron) and Q 2 from a flat distribution, and then calcu- 
lated all remaining momenta and angles constrained by 
quasi-elastic kinematics. If the electron and the proton 
fell in the CLAS acceptance we recorded the kinematics 
of the event in a form of an Ntuple Q and we weighted 
it with the corresponding cross section, differential in 
p n , n , 4> n , Q 2 and (f> e . The events were then binned 
identically to the experimental data using the same cuts. 
No normalization factors between theoretical and exper- 
imental data were used. 

This model is an extension of earlier diagrammatic 
methods [ToL 11 1 to JLab kinematics. It incorporates 
four amplitudes: the Plane Wave Impulse Approxima- 



tion (PWIA), meson exchange currents (MECs), high en- 
ergy diffractive nucleon-nucleon elastic scattering (FSIs) 
and intermediate A-nucleon rescattering (AN). Deuteron 
wave functions derived from both the Paris [l2| and 
the Argonne V18 [l]| potentials were used. The elec- 
tron couples to the nucleons through a fully relativis- 
tic, on-shell nucleon current. The dipole parameteriza- 
tion was chosen for the magnetic form factors of the nu- 
cleon. The latest JLab data 14J were used for the proton 
electric form factor, while the Galster [l5| parameteriza- 
tion was selected for the neutron electric form factor. 
The parameters of the NN amplitude are the same as in 
Ref. [8J], and are fixed by the elastic scattering cross sec- 
tion. The 7r and p exchanges are taken into account in 
the MEC and AN formation amplitudes, as described in 
Ref. llj. The electromagnetic N — > A transition form 
factor F NA {Q 2 ) = (1 - Q 2 /9)/(l + Q 2 /0.Tj 2 is driven by 
the world data (MAID parameterizatio n [la] ) and specif- 
ically by the highest Q 2 measurement [17] in Hall C at 
JLab. The most recent data set [l8[ from CLAS is lower 
by as much as 10% for Q 2 < 3 GeV 2 but is similar for 
Q 2 > 3 GeV 2 . 

The calculated cross sections are shown in Figs.[l][2j[3l 
and [4] Systematic uncertainties in the theoretical cross 
sections come from the on-shell approximation for the 
electron- nucleon current (~5%), the parameterization of 
the NN elastic scattering amplitude (~10%), and the pa- 
rameterization of F/va (~11%). Thus, the systematic 
uncertainties in the theoretical predictions are ~ 15% for 
the full calculations. Since the MEC amplitude in our Q 2 
range is very small, the corresponding uncertainty can be 
neglected. 

Figs. Q] and [2] show the distributions in recoil neutron 
momentum integrated over the angular range 20° < 9n< 
160°, where acceptance corrections are well defined [6J. 
The experimental p n distribution drops by three orders 
of magnitude over the range 0-2 GeV/c similar to the full 
theoretical calculations. For p n < 800 MeV/c, however, 
the data and calculations agree better than for higher 
p n - Below p n — 250 MeV/c, quasi-elastic scattering of 
electrons on protons (the PWIA channel) exhausts the 
cross section. Neutron-proton FSI dominates for 250 < 
p n < 750 MeV/c, while intermediate A production is 
prominent for p n > 750 MeV/c, bringing the model close 
to the data. Both Paris [13] and Argonne V18 [l]| wave 
functions show similar results for p n < 1 GeV/c, whereas 
above 1 GeV/c the two wave functions differ strongly, 
and lead to very different PWIA contributions. However, 
the AN channel overwhelms the cross section here: low 
momentum components of the wave function feed these 
higher values ofp n , and the sensitivity of the cross section 
to the high momentum components of the wave function 
is lost. Nevertheless, the theory agrees well with the data. 

Although the theory describes the neutron momentum 
distributions well, the log scale makes a close comparison 
difficult. The remaining differences between theory and 
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FIG. 1: Color online. The recoil neutron momentum dis- 
tribution for (a) Q 2 = 2 ± 0.25 GeV 2 and (b) Q 2 = 3 ± 
0.5 GeV 2 . Dashed, dash-dotted and solid curves are calcu- 
lations with the Paris potential for PWIA, PWIA+FSI and 
PWIA+FSI+MEC+NA, respectively. Dotted (red) curves 
are calculations with the AV18 potential. 



experiment are best seen quantitatively in the linear plots 
of angular distributions for various regions in p n below 
600 MeV/c. The detailed theoretical investigations in 
Refs. [13, [13] have shown that there are specific features 
of recoil neutron angular distributions for different ranges 
in p n . For p n < 0.1 GeV/c the angular distributions are 
expected to be insensitive to FSIs, for p n ~ 0.4 — 0.5 
GeV/c FSIs should dominate, and for 0.2 < p n < 0.3 
GeV/c the interference between PWIA and FSI ampli- 
tudes should contribute. Until now, this characteristic 
behavior of the recoil neutron angular distributions has 
not been checked experimentally. 

Figs. [3] and [4] show neutron angular distributions for 
three ranges of p n at Q 2 — 2, 3, 4 and 5 GeV 2 . Each panel 
clearly shows the evolution of the interaction effects with 
p n and for a fixed value of Q 2 . 

In the highest momentum range (0.4 < p n < 0.6 
GeV/c) the angular distributions exhibit a large peak in 
the vicinity of 6 n = 70°. This effect comes from neutron- 
proton rescattering, and corresponds to the on-shell prop- 
agation of the struck nucleoli. It is maximal when the 
kinematics allow for rescattering on a nucleoli almost at 
rest [13], which happens when x = Q 2 jlMv = 1 {v is the 
energy of the virtual photon, and M is the nucleon mass). 
The following physical picture emerges. The electron 



FIG. 2: Color online. The same as Fig(T] for (a) Q 2 = 4 ± 
0.5 GeV 2 and (b) Q 2 = 5 ± 0.5 GeV 2 . 



scatters primarily from a proton almost at rest. Since 
the total energy is larger than the sum of the masses of 
the two nucleons, the struck proton can propagate on- 
shell and rescatter off the neutron which is also nearly 
at rest. In the lab frame, the soft neutron recoils at 90° 
with respect to the fast forward proton. Two-body kine- 
matics places the rescattering peak at about 8 n =70° for 
our kinematics. In the classical Glauber approximation, 
the nucleon propagator is linearized and recoil effects are 
neglected, and herefore, the rescattering peak stays at 
n =90° [U E2]. This has been fixed in the General- 
ized Eikonal Approximation (GEA) [13] which takes into 
account higher order recoil terms in the nucleon propa- 
gator. In the diagrammatic approach the full kinematics 
are taken into account from the beginning [3, [13] ■ The 
shape of the angular distribution reflects the momentum 
distribution of the proton in deuterium. 

A A resonance produced on a nucleon at rest at x = 
[1 + (Ml - M 2 )/Q 2 ] -1 < 1, can propagate on-shell and 
rescatters from the second nucleon also at rest ||. This 
contribution shifts the rescattering peak toward larger 
angles, and brings the theory into better agreement with 
experiment. It also decreases faster with Q 2 , consistent 
with the steeper variation of the N — * A transition elec- 
tromagnetic form factor as compared to the dipole pa- 
rameterization of the nucleon form factors. The excess 
theoretical cross section at Q 2 = 2 GeV 2 is a reflection 
of our linear fit to the ratio of N — > A and dipole form 
factors. A better fit to the latest data [13] from CLAS 
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FIG. 3: Color online. The recoil neutron angular distribution 
for (a) Q 2 = 2 ±0.25 GeV 2 , 400 < p n < 600 MeV/c; (b) Q 2 = 
3 ±0.5 GeV 2 400 < p n < 600 MeV/c; (c) Q 2 = 2 ±0.25 GeV 2 , 
200 < p n < 300 MeV/c; and (d) Q 2 = 3 ± 0.5 GeV 2 , 200 < 
p n < 300 MeV/c. The data for p n < 100 MeV/c are plotted 
in the bottom part of (c) and (d) and scaled by 0.035. The 
curves have the same meaning as in Fig. [1] 



FIG. 4: Color online. The same as in Fig [3] but for Q 2 = 
4 ± 0.5 GeV 2 ((a) and (c)) and Q 2 = 5 ± 0.5 GeV 2 ((b) and 
(d))- 



leads to a reduction of the peak by ~ 15 % for Q 2 < 3 
GeV 2 , in better agreement with experiment. 

In the intermediate momentum range (0.2 < p n < 0.3 
GeV/c), FSIs suppress the quasi-elastic contribution in 
the vicinity of x = 1. Here the relative kinetic energy be- 
tween the outgoing proton and neutron T ~ Q 2 /2M lies 
between 1 and 3 GeV. The nucleon-nucleon scattering 
amplitude is almost purely absorptive and the FSI am- 
plitude interferes destructively with the quasi-free ampli- 
tude. This induces a loss of flux for fast protons. 

In the lowest momentum range (p n < 0.1 GeV/c) 
rescattering effects are small, and the experimental and 
theoretical angular distributions are similarly flat. The 
magnitude of the experimental cross sections is well re- 
produced at low Q 2 , but the theory slightly exceeds the 
data at larger Q 2 . This effect has already been ob- 
served in the study of 3 He(e,e'p) 2 H at low recoil mo- 
mentum [23} , and is not yet well understood. 

In summary, this benchmark experiment demonstrates 
that the mechanisms of the exclusive 2 H(e, e'p)n reac- 
tion are well understood for 1.75 < Q 2 < 5.5 GeV 2 . 
Theoretical and experimental cross sections agree within 
20%, consistent with the systematic uncertainties («15% 



for theory and «10% for experiment). Proton-neutron 
rescattering (FSIs) and A production dominate over a 
large part of the phase space, except at backward an- 
gles (8 n > 110°) or very low recoil momenta (p n < 
100 MeV/c), where the distributions directly reflect the 
deuteron wave function. A good understanding of the 
mechanisms of the exclusive 2 H(e, e'p)n reaction in our 
Q 2 region opens an opportunity to investigate the short 
distance properties of nucleons in deuterium, which will 
be discussed in our future publications. 
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